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Abstract 

We study a generalised model of population growth in which the state variable is 
population growth rate instead of population size. Stochastic parametric perturba- 
tions, modelling phenotypic variability, lead to a Langevin system with two sources 
of multiplicative noise. The stationary probability distributions have two charac- 
teristic power-law scales. Numerical simulations show that noise suppresses the ex- 
plosion of the growth rate which occurs in the deterministic counterpart. Instead, 
in different parameter regimes populations will grow with "anomalous" stochastic 
rates and (i) stabilise at "random carrying capacities" , or (ii) go extinct in random 
times. Using logistic fits to reconstruct the simulated data, we find that even highly 
significant estimations do not recover or reflect information about the determin- 
istic part of the process. Therefore, the logistic interpretation is not biologically 
meaningful. These results have implications for distinct model-aided calculations 
in biological situations because these kinds of estimations could lead to spurious 
conclusions. 
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Population dynamics are frequently modelled with simple equations that mimic 
some aspects of replicating biological entities, such as division (in cells), fis- 
sion (in modular organisms) or reproduction (in eukaryotes), competition, 
and population-size limiting (saturation). These and other properties are rep- 
resented by various models. Frequently the validity of these models is a matter 
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of statistical goodness of fit with a specific data set. However, tliese biolog- 
ical properties are not entirely of intrinsic nature to the individuals, or to 
the populations themselves, but rather emerging ecological properties, i.e. the 
interaction between "individuals" and "environment" . The models of popula- 
tion growth simplify (whenever it is possible) the potential complexity of a 
detailed ecological description into simple equations. 

Recently we showed that a variety of biological growth models can be unified 
using a phase-space decomposition using two dynamical variables, population 
size X and growth rate r [2] (analogous to a particle' s position and momentum, 
respectively) : 

X = xr , (la) 
r = r{6r — p). (lb) 



The constant p is the Malthusian parameter, and 6 is the intraspecific inter- 
action coefficient. By varying these two parameters it is possible to reproduce 
exactly a wide family of growth laws including exponentials, logistics [il. E3|. 
Gompertzian 0], Potential [1, 01, as well as allometric growth laws like Von 
Bertalanffy's [JIBl and West's equations, among others [l|. 



Although the deterministic behaviours can be associated with distinct biolog- 
ical scenarios jl|, populations are often influenced by some source of noise. 
Biologically, random "forces" are often related to environmental and demo- 
graphic fluctuations, as well as to intrinsic complex effects like genetic variabil- 
ity, mating, and segregation. These environmental and genetic factors can be 
thought to be the determinants of the parameters that describe the growth of 
a population. This conjecture means that they are (complicated) measures of 
phenotypic expression. Thus to a flrst approximation we can model the effect 
of phenotypic variation as noise over these parameters. 

We can consider two sources of noise rjp (j) and rjg (t) affecting respectively 
the parameters 6 and p. The rate equation (lb) linearly perturbed with p — >■ 
P + Vp (i) and + Tje {t) results in 

r = r{er - p) + r'^rjg ^t) + rrjp ^t) , (2) 

where the rji (j) have the usual properties of white noise: 

iVi (t)) = 

{Vi {t)Vi is)) = S{t - s)e^i (3) 
iVe it)Vp it)) =7 
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with (. . .) denoting expectations, and i = p,9. Here, are the intensities of 
the noise sources, 7 is the correlation between the two noise sources, and 6(t) 
is Dirac's delta function. 



The resulting system is a Langevin equation where p is the drift term and 
can be thought as the force of an external field j9|. Multiplicative noise, 
often represents fluctuating barriers or processes of anomalous diffusion (i.e. 
diffusion where the probability of long steps is higher then in the normal case) 
[1^ 11, l35 0,E3- Also, multiplicative noise is a process that retains memory 



(i.e. is non-Markovian) , and has been investigated in the context of population 
growth and extinctions H, 3, [l^ ■ 



The stochastic differential equation (SDE, Eq. 2) remains uncoupled from the 
size X. This gives an operational advantage since the analyses of the SDE can 
be made in terms of r as a 1-dimensional system that is relatively simple to 
handle. 

To study the effects of multiplicative noise, and make precise the meaning of 
"anomalous growth" in populations (in analogy to anomalous diffusion), first 
consider the probability distribution for the rates. In the Ito interpretation of 
noise the probability is given by the related Fokker Planck equation (FPE): 
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r{er - p)V^r,t) +^drr ie;-2epee-fr + eyyV(^r,t) ■ (4) 



Setting the time derivative equal to zero makes it possible to calculate the 
potential solution of the FPE on the stationary regime P(r) (i.e. equilibrium 
solution), which gives 

P(r) := Af (r-2)^«+^ (eg - eo^fr + rY'~' x 
X exp[2(ec6'o-7cPo)tan ^ (e^r - 7^)] , 



where J\f is the integration constant, and 

Po = p/ep, eo = e/el, eo = ep/ee 

ec = 6o^(l-7^)-V^7c = 7(l-7^)-^/^■ 



Fig. 1 shows that there are three distinct kinds of stationary distributions for 
the rates. The first thing that we note, is that the correlation 7 modulates 
the transition from one distribution to another. Thus, for simplicity for the 
further analyses we proceed setting 7 = 0. 
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The first distribution (Fig. lA) is monotonous decreasing, and the other two 
(Figs. IB-C) have an analytic maximum at 



'9n 



- 4{po + 1) 



The condition for having an analytic maximum is 

I > elipo + 1) . (6) 



Fig. 2 outlines the regions where the inequality (6) holds. The parabolic curve 
(i.e. the boundary given by the equality in Eq. 6) divides the parameter space 
into three regions with different stationary regimes. The first, when < 6q/8 < 
eo(po + l), corresponds to the space under the parabola (Fig. 2); in this case the 
probability density is accumulated at r = 0. The second region, characterised 
by < el{po + 1) < Oq/8, corresponds to the space over or under the parabola 
region (Fig. 2). The third region, defined by eQ(po + l) < < 0^/8, is the space 
on the left of the parabola. In the two last regions, the probability mass of rates 
is distributed along the axis, indicating that the growth rate can be persistent 
(i.e. non-zero). In the following of the paper, we will show that each of these 
regions have distinct qualitative solutions in which the deterministic nature of 
the process is "forgotten", but the resulting dynamics of the population size 
look like exponential or logistic dynamics. We will demonstrate however, that 
these two forms are entirely product of noise, hence fitting these models to 
the realizations -although statistically significant- are spurious. 



For particular cases of Eq. (2) the stationary distribution has been calculated 
before. When 6 and its nois e m term are absent, the equation recovers the 



geometric Brownian motion [18[, whose stationary distributions were shown 
to have power- law tails 0. I~h,s represe„tat,o„: the gxowth corresponds to 
a Gompertzian growth. A power-law-tailed distribution is also found for the 
stationary distribution of an equation where p 7^ but which is not perturbed 
[lol . Also, a logistic case 6 = 1 was analysed by Morita and Makino [2^ using 
perturbation techniques for a time dependent solution. In log-log scale the 
distribution (5) is kinked near Tc = eoexp (7r6'o/2eo(l — Po)), with a right-tail 
decreasing in a power law fashion log P(r) ~ — 41ogr (insets in Fig. 1). This is 
a result that can be derived directly from the particular case studied by Gora 
[1^ , because the right tail of the distribution is independent of the parameters. 
Moreover, there is also a power law behaviour for small values of r which is 
given by log P(r) ~ — 2(po + l)logr. These power law tails lead to Tsallis 
statistics [lO|, |21|. Some relationships between exponents have been derived 
for a system related (but not equivalent) to ours |22i |. 
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The fast decrease of the right tail has an important consequence, which is the 
boundedness of the process. In other words, it means that the fluctuations 
remain finite. The simulations of Fig. 3A show that when the probability 
density is accumulated at r = the rates will stochastically reach zero and 
stay there forever. Whenever this happens, the population freezes at its - 
random- current size. 



Sibly et al. performed an analysis where they fitted more than a thousand 
population time series to the 6'-logistic model. Their analysis was based on 
the size-dependent per capita growth rate r{x). The examples they presented, 
show comparable patterns to the realizations obtained from our model (Fig. 
3A). However, as we can see in this figure, the stochastic trajectories are not- 
centredd on the deterministic trajectories, as it is common for multiplicative 



perturbations. Therefore, the interpretation of the estimations in Ref. |24[ dif- 
fer from the deterministic path, at least in the light of our model. We will 
return to this discussion later in the article. 

Other simulations, using a processes having stationary distribution with max- 
ima, are shown in figure 3B-C. In this cases the rate does not explode in the 
time-window, even when these realizations (the deterministic and the stochas- 
tic) have the same initial conditions which would lead to explosions in the 
absence of perturbations. 

Recently, Mao et al. 2^ demonstrated that the deterministic explosions of 
"positive" logistic equations (e.g. of the form x = ax{l + bx)) can be controlled 
with certain types of multiplicative noise sources. When these fluctuations 
are present, populations will not diverge in finite time, although their purely 
deterministic analogue does. The rate-representation introduced in this paper 
is also of quadratic form, thus the results of Mao et al. apply to Eq. 
(2). However, the biological interpretations change, because explosions are 
suppressed in the rate rather than in population size. 

As indicated by the distribution of the rates, probability is accumulated near 
the maximum, thus the rates will be non- vanishing, jumping from very slow to 
high (but finite) values, making the size of the population increase in bursts, 
reconstructing a devil's staircase pattern (a staircase where all the steps are 
of different size and height). Also, because the rate never reaches zero, the 
population grows unlimited. 

The same distributions of Fig. 1 appear for negative rates. The course of 
the population is the opposite, i.e. decreasing, although the distribution is the 
same (in absolute value): (i)if the distributions have an analytic maximum, the 
rates will remain finite and fluctuating, meaning that population will decrease 
erratically but monotonously and therefore populations will become extinct in 
random times; (ii) if the distributions do not have an analytic maximum, then 
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r reaches zero stochastically (Fig. 3A), and then the populations will stabilise, 
again at a "random carrying capacity" . 

The distribution (5) is not normalisable whenever po + 1 > 0, because it 
diverges when r 0. This means that r = is an exit barrier, and hence once 
the rate reaches zero it will stay there. This limit is the same if taken from the 
left, thus the rate cannot either jump to a negative value once it reaches zero. 
For instance the rates maintain their sign or become null, but never change 
sign. The meaning is that when the rates are stationary, an initially growing 
populations will continue to grow, or at most, cease growing but they will not 
suddenly shrink. Therefore, the converse is also true: populations that started 
shrinking, will not suddenly change its course and grow. They will continue 
to shrink until extinction, or reach a stable value. 

The dynamics results in distinct realizations that can give drastically different 
solutions, when compared, for example, to the deterministic solution. Fig. 
4A show that the equilibrium value of the populations can be very different 
from the deterministic carrying capacity. Thus the observed equilibrium value 
of the populations is no longer determined by the initial conditions, as in the 
deterministic case jfj. Actually, the carrying capacity is now a random variable. 
For example, in Fig. 4A the size equation is solved for several realizations of 
the process (2). For the naive eye, the distinct realizations could be seen as 
distinct "noisy logistics" with different carrying capacities. Comparing the 
realizations to a logistic equation gives highly significant fits, even when the 
data come from a common process having the same values of the parameters. 

In order to determine if we can recover deterministic information of the pro- 
cesses, we performed simulations of 250 randomly selected values of p, 9 and 
for each we performed 30 realizations. To every growth curve we least-squares- 
fitted a logistic model, and calculated its parameters p and x^o- Fig. 4B shows 
a scatter plot of the estimated vs. the deterministic values of the Malthusian 
parameter, showing a poor relationship. Other correlations are presented in 
Table 1. These results show that the reconstructions are totally spurious since 
they do not reflect any information of the generating process. But because the 
dynamics resemble a logistic realization, accepting a null hypothesis that the 
biological phenomena determining growth are of logistic nature is true, sta- 
tistically speaking. However, our calculations show that stochastic processes 
can account for the same qualitative and quantitative description. Therefore, 
simple analysis like least squares fits are not enough to confirm the logistic 
hypothesis. 

There is an analogous effect for the case when the rates are persistent. Once the 
stochastic rates are in stationarity, the resulting population dynamics resem- 
bles exponential growth. In the deterministic exponential growth, the initial 
condition of the rate determines the growth parameter 0]. However, under our 
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scheme, the growth process is Markovian, thus the initial conditions of (x, r) 
do not affect the stationary distribution. As a consequence, the expected or 
averaged rates are spurious estimators of an exponential dynamic (Fig. 4C). 
A similar problem was described by Renshaw [26] when demographic stochas- 
ticity is present in an exponentially growing population. 

At this point it is necessary to make a distinction between the outcomes of 
noise sources coming from demographic or phenotypic stochasticity. The first 
has been studied and experimentally supported [la ■ This kind of stochas- 
ticity is such that randomness affects the population through events of acci- 
dental mortality or occasional migrations (and is analogous to energy input 
coming from a heat bath.) In these cases, the populations would fluctuate, for 
example, close to carrying capacities, and thus information for the determin- 
istic part of the dynamics can be extracted by averaging. The second type, i.e. 
parameter stochasticity, is more related to fluctuations in phenotypes, which 
results from the "superposition" of genetic and environmental processes. How- 
ever, from the perspective of our model, where carrying capacities are not an 
intrinsic property of the environment, this averaging might not make biologi- 
cal sense. As we said, an average of the stochastic trajectory does not recover 
the deterministic path, like in Figs. 4A, 4A. Of course, populations might 
still be subject to demographic stochasticity, and therefore show fluctuations 
around a stable size. In this case, we would be presented with an additional 
noise source r]M, more related to the measuring techniques, perturbing the size 
equation as: 

X = xr + rjM , (7) 



that gives the fluctuating pattern over the stable size. (This is a problem known 
as filtering: when the measuring procedure has additional noise sources, not 
taking them into account in the estimations, may bias the interpretation of 
the underlying process [3].) Considering this source of fluctuations is more 
related to time series estimation than to the biological aspects of our model. 
Therefore we shall not discuss it in further detail, and we defere the reader to 
Ref. ji^ for a method to overcome these kind of problems. 

To summarise, we have presented an analysis of a novel population growth 
model that is based in fluctuations in the per-capita growth rate, rather than 
in the growth variable. The result, is that the rate always remains finite, ei- 
ther because rate explosions are suppressed, or because rate is damped to zero 
stochastically. As a consequence, and depending on the relationship between 
the deterministic parameters and noise, the model reproduce patterns that 
resemble exponential and logistic (sigmoid) growth. It is important to notice 
that these behaviours are irrespective on how the deterministic population 
would grow. These forms are determined by the fluctuations and not from 
the biological processes of birth and death, at least not in the conventional 
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interpret at ionn and description. When ^q/8 > eQ(po + 1) then the resulting 
population grows anomalously, but with bounded fluctuations, and resem- 
bles an exponential growth. When 6*0/8 < el{po + 1) then the populations 
grow toward saturation. However, this result challenges the idea of a carrying 
capacity, that is supposed to describe self regulatory processes and an intrin- 
sic property of the environment. Here, it is an emerging property from the 
fluctuations. In both cases, and more critically in the second (the logistic), 
statistical fits to the realizations are highly significant. But since the effects 
of randomness override the deterministic forces of the system, making these 
statistical estimations becomes unreliable. In the context of our formulation, 
the question about rate estimations looses its sense, because forecasting us- 
ing the classic deterministic models proves useless. Thus fluctuation analysis 
might prove more informative about the stochastic driving forces. In this way, 
estimations and forecasting can give other statistical solutions to classical and 
new problems, using our different perspective, that is, when populations are 
subject to phenotypic stochastic variability. 
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Table 1 

The correlation coefficient between the parameters and estimators for the logistic 
estimations are very low, indicating that the fits have little predictive power with re- 
spect to the deterministic part of the process, even thought they could be considered 
virtually perfect. 





P 


e 


P 


p 


0.532 


0.301 






-0.0712 


-0.0290 


-0.0119 



10 



Figure Captions. 



Fig. 1. Potential solutions for the Stationary Fokker Planck equation. (A) When 
< ^q/8 < eQ(/9o + 1), the distribution shows an exit barrier at r = 0; in this 
example p = 1, Cp = 0.5, 9 = 1, ee = 0.1. (B) When < eg(po + 1) < 6*^/8, 
besides the exit barrier, the distribution also shows an analytic maximum; in this 
example p = 1, €p = 0.4, 9 = 4, €0 = 1. (C) When eg(po + 1) < < 9^/8, the exit 
barrier disappearss, and the distribution is unimodal with an analytic maximum; the 
parameters are p = —3, €p = 0.5, 9 = 2, €$ = 1.5. The insets plot the distributions 
in log-log scale, showing that there are two characteristic scales. The left tail scales 
with an exponent of —2(1 + po) whereas the right tail scales with an exponent of 
—4. The inflection points are close to the maximum. 

Fig. 2. The parameter space (p, 9) consists of four quadrants, corresponding to 
their sign combinations. The dotted lines indicate distinct deterministic growth 
functions known in the literature: Potential (P) p = 0; logistic (L) p > 0,9 = I; 
Gompertzian (G) p > 0, 6* = 0; West (W) p > 0,9 = -1/4; von Bertalanffy (VB) 
p > 0, ^ = —1/3; Exponential (E) p = 0, ^ = 0. The solid curve represents the 
noise-transition points between the three distinct regimes of the distributions of the 
rate: (a) inside the parabolic region < 6*0/8 < €q{po + 1); (b) above or below the 
parabolic region < e^ipo + 1) < ^o/^! (c) at the left of the parabolic region 
elipo + 1) < < 9i/8. 

Fig. 3. Realization of noisy rate dynamics. The dotted lines show deterministic 
solutions, while the continuous bold lines show the stochastic realizations. It can 
be seen that when the deterministic rates decrease to zero, the stochastic dynamics 
will also decrease to zero. Also, when deterministic rates explode the stochastic 
dynamics remain finite. These realizations correspond respectively to the regimes 
and parameters of Fig. 1 
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Fig. 4. (A) Integrations for population size for distinct realizations of the same pro- 
cess when 0q/8 < eo(/3o + 1) using the parameters p = 1, Cp = 0.5, = l,ee = 0.5. 
The bold line shows the deterministic dynamics and the thin lines are realization 
for the population size. The dotted lines are logistic estimations. In this cases, the 
estimated Malthusian parameters range between p G (0.8, 1.5). Carrying capacities 
range from x^^i G (1, 10^), and the deterministic carrying capacity is Xoo = 10.0. All 
of the estimations have a regression coefficient > 0.995 with p-values less than 
10"'^. (B) Correlation between the estimated and generating Malthusian param- 
eters, from 7500 simulation spanningg 250 distinct pairs of uniformly distributed 
values of p G [0,2], 9 = ±eo^/2{^^Tl), using Cp = 1.0, eg = 0.1. The radii of the 
circles are 10~^ log(S'£^) («S'-E' = standard error). The continuous line is the linear 
trend, which gives p = 0.389591 + 0.371044p with i?^ = 0.2873 (p< 2.2 lO^^^). 
The dotted line is a linear trend weighted with the inverse of the standard error 
of each estimation: p = 0.91300 + 0.22528p with i?^ = 0.3393 (p< 2.2 10"^^). 
Comparing these two estimations we see that even in the best case (the weighted 
regression) the predictive power is poor. (C) Integration for population sizes for 
distinct realizations of the process when 6^/8 > e^po + 1) using the parameters 
p = —3, Cp = 0.3, 9 = 2,eg = 1.5). The bold line represents the deterministic dy- 
namics, and the thin lines the realizations for population size. The dotted lines are 
estimations for the exponential growth. The estimated values for the exponential 
growth parameter are in the range of (21.8, 26.0). All the estimations have a regres- 
sion coefficient > 0.998 with p-values less than 10~^^. The graph is in semi-log 
scale. 
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